********************************************************************************
/*
							Master do-file for
			
					Norms, Enforcement and Tax Evasion

		     by Tim Besley, Anders Jensen and Torsten Persson

			 
This do-file creates all the figures and tables contained in the main text
and the appendix for the article that is forthcoming in Review of Economics and
Statistics

*/
********************************************************************************

clear
clear matrix
* Set directory to location where replication package is downloaded
cd "C:\Users\ajensen\Dropbox (CID)\Research\projects\T T project\Submission_April 2020\ReStat submission\RR_R1\Resubmitted\Replication Package_Besley et al_2021" 

use "Data sets\dwellings per band_1993 2009.dta"
sort borough year
egen group=group(borough year)
egen unique=count(year), by(group)
tab unique
drop if unique>1 & unique!=.
drop unique group
sort borough year

merge 1:1 borough year using "Data sets\public finance and vote data.dta"
drop if _merge==1

drop _merge

merge 1:1 borough year using "Data sets\tax liab diff bands.dta" 
drop if _merge==2

drop _merge
merge 1:1 borough year using "Data sets\average earnings per capita.dta"

gen lg_dwelincome=log(inc)
replace lgincome=lg_dwelincome if year>=2002 & _merge==3
drop if _merge==2
drop _merge
drop inc

merge 1:1 borough year using "Data sets\popl.dta"

gen logpop=log(popl)
replace lg_pop=logpop if year>=1999 & _merge==3
drop logpop
drop _merge

egen group=group(bor_id year)
egen unique=count(year), by(group)
tab unique
drop if unique>1 & unique!=.

drop unique group

xtset bor_id year

local dwel_stock "dwellings_f dwellings_g dwellings_h"
foreach var of local dwel_stock{
ipolate `var' year if year<=1996, gen(`var'_interpol)
replace `var'=`var'_interpol if `var'==.
}
gen dwellings_total=dwellings_a+dwellings_b+dwellings_c+dwellings_d+dwellings_e+dwellings_f+dwellings_g+dwellings_h if dwellings_a!=. & dwellings_b!=. & dwellings_c!=. & dwellings_d!=. & dwellings_e!=. & dwellings_f!=. & dwellings_g!=. & dwellings_h!=.

local dwelling_band "dwellings_a dwellings_b dwellings_c dwellings_d dwellings_e dwellings_f dwellings_g dwellings_h"
foreach var of local dwelling_band{
gen sharetotl_`var'=`var'/dwellings_total
}

sort year

gen h_to_a=sharetotl_dwellings_h/sharetotl_dwellings_a
egen htoa1993=mean(h_to_a) if year==1993, by(bor_id)
egen h_to_a_1993=mean(htoa1993), by(bor_id)

gen sharehh_renters_locauth=totalhhrentedlocauth/number_dwellings89
egen sharehhrentlocauth=mean(sharehh_renters_locauth), by(bor_id)

gen sharehh_renters1=totalhhrentedprivatefurnished/number_dwellings89
gen sharehh_renters2=totalhhrentedprivateunfur/number_dwellings89
egen sharehhrent1=mean(sharehh_renters1), by(bor_id)
egen sharehhrent2=mean(sharehh_renters2), by(bor_id)
gen sharepriv=sharehhrent1+sharehhrent2

sort year
********************************************************************************
*********Figure 1 A **
egen meanyear_ct_evasion=mean(ct_evasion_allyears), by(year)
twoway connected  meanyear_ct_evasion year, mcolor(blue) msymbol(circle) lcolor(blue) lwidth(medthick) msize(small) xline(1990 1992) xlabel(1980(5)2005) xtitle("Year") ytitle("Percent") text(10 1998 "1990-1992: Poll Tax base") graphregion(color(white))
graph save "Figures and Tables\Fig1A.gph", replace
********************************************************************************

********************************************************************************
************** Figure 1B **
kdensity ct_evasion_allyears if year<1989 & ct_evasion_allyears<=30, lpattern(solid) lcolor(black) plot(kdensity ct_evasion_allyears if year>=1990 & year<=1992 & ct_evasion_allyears<=30, lpattern(shortdash) lcolor(red)  || kdensity ct_evasion_allyears if year>=1993 & year<=1994 & ct_evasion_allyears<=30, lpattern(solid) lcolor(blue) || kdensity ct_evasion_allyears if year>=1995 & ct_evasion_allyears<=30, lpattern(longdash) lcolor(dkgreen)) legend(label(1 "1980-89") label(2 "1990-92") label(3 "1993-94") label(4 "1995-09")) legend(cols(4)) xtitle("Tax Evasion %") ytitle("Marginal Density") title("") graphregion(color(white))
graph save "Figures and Tables\Fig1B.gph", replace
********************************************************************************

********************************************************************************
egen mean_ptevasion=mean(ct_evasion_allyears) if year>=1990 & year<=1992, by(bor_id)
egen mean_ptevasion_council=mean(mean_ptevasion), by(bor_id)
egen median_ptevasion=median(mean_ptevasion)

egen mean_drevasion=mean(ct_evasion_allyears) if year<1990
egen mean_drevasion_all=mean(mean_drevasion)
egen mean_drevasioncouncil=mean(ct_evasion_allyears) if year<1990, by(bor_id)
egen mean_drevasion_council=mean(mean_drevasion_all), by(bor_id)
gen mean_ptdrevasion_council=mean_ptevasion_council-mean_drevasion_all
egen median_ptdrevasion=median(mean_ptdrevasion_council)

egen mean_drevasionshort=mean(ct_evasion_allyears) if year<1990 & year>=1987
egen mean_drevasionshort_all=mean(mean_drevasionshort)
gen mean_ptdrevasionshort_council=mean_ptevasion_council-mean_drevasionshort_all

gen high_evasion=.
replace high_evasion=1 if mean_ptdrevasion_council>median_ptdrevasion & mean_ptdrevasion_council!=.
replace high_evasion=0 if mean_ptdrevasion_council<=median_ptdrevasion & mean_ptdrevasion_council!=.

********************************************************************************
************** Figure 2A *******************************************************
sort year
egen mean_ctevasion_lowevasion=mean(ct_evasion_allyears) if high_evasion==0, by(year)
egen mean_ctevasion_highevasion=mean(ct_evasion_allyears) if high_evasion==1, by(year)
twoway (connected mean_ctevasion_lowevasion year, lcolor(blue) lwidth(medthick) msymbol(circle) msize(small) mcolor(blue)) (connected mean_ctevasion_highevasion year, msymbol(circle) msize(small) mcolor(red) lcolor(red) lwidth(medthick)), xline(1990 1992) legend(label(1 "PT Evasion<Median") label(2 "PT Evasion>=Median")) xtitle("Year") ytitle("Percent") text(16 1998 "1990-92: Poll Tax base") graphregion(color(white))
graph save "Figures and Tables\Fig2A.gph", replace
********************************************************************************

**************************************
*** Conditional poll tax evasion *****
xi: reg ct_evasion_allyears ct_taxliability_dwellin i.region year_dum11 year_dum13 if year>=1990 & year<=1992, robust
xi: reg ct_evasion_allyears ct_taxliability_dwellin lgincome i.region year_dum11 year_dum13 if year>=1990 & year<=1992, robust
xi: reg ct_evasion_allyears ct_taxliability_dwellin lgincome lg_pop i.region year_dum11 year_dum13 if year>=1990 & year<=1992, robust
xi: reg ct_evasion_allyears ct_taxliability_dwellin lgincome lg_pop con_share conctl lab_share labctl i.region year_dum11 year_dum13 if year>=1990 & year<=1992, robust
xi: reg ct_evasion_allyears ct_taxliability_dwellin lgincome lg_pop con_share conctl lab_share labctl h_to_a_1993 sharepriv sharehhrentlocauth i.region year_dum11 year_dum13 if year>=1990 & year<=1992, robust
predict res if year>=1990 & year<=1992, residuals
egen mean_ptresidual=mean(res), by(bor_id)
gen mean_ptresdr_council=mean_ptresidual-mean_drevasion_all
egen median_ptresdr=median(mean_ptresdr_council)
egen median_ptres=median(mean_ptresidual)
gen mean_ptresdrshort_council=mean_ptresidual-mean_drevasionshort_all

gen high_res=.
replace high_res=1 if mean_ptresdr_council>median_ptresdr & mean_ptresdr_council!=.
replace high_res=0 if mean_ptresdr_council<=median_ptresdr & mean_ptresdr_council!=.
**************************************

************************Fig 2B**************************************************
egen mean_ctevasion_lowres=mean(ct_evasion_allyears) if high_res==0, by(year)
egen mean_ctevasion_highres=mean(ct_evasion_allyears) if high_res==1, by(year)
sort year
twoway (connected mean_ctevasion_lowres year, msymbol(circle) msize(small) mcolor(blue) lcolor(blue) lwidth(medthick)) (connected mean_ctevasion_highres year, msymbol(circle) mcolor(red) msize(small) lcolor(red) lwidth(medthick)), xline(1990 1992) legend(label(1 "PT Conditional Evasion<Median") label(2 "PT Conditional Evasion>=Median")) xtitle("Year") ytitle("Percent") text(18 1998 "1990-1992: Poll Tax base") graphregion(color(white))
graph save "Figures and Tables\Fig2B.gph", replace
********************************************************************************

gen hres_1993=high_res*year_dum14
gen hres_1994=high_res*year_dum15
gen hres_1995=high_res*year_dum16
gen hres_1996=high_res*year_dum17
gen hres_1997=high_res*year_dum18
gen hres_1998=high_res*year_dum19
gen hres_1999=high_res*year_dum20
gen hres_2000=high_res*year_dum21
gen hres_2001=high_res*year_dum22
gen hres_2002=high_res*year_dum23
gen hres_2003=high_res*year_dum24
gen hres_2004=high_res*year_dum25
gen hres_2005=high_res*year_dum26
gen hres_2006=high_res*year_dum27
gen hres_2007=high_res*year_dum28
gen hres_2008=high_res*year_dum29
gen hres_2009=high_res*year_dum30

gen hres_1980=high_res*year_dum1
gen hres_1981=high_res*year_dum2
gen hres_1982=high_res*year_dum3
gen hres_1983=high_res*year_dum4
gen hres_1984=high_res*year_dum5
gen hres_1985=high_res*year_dum6
gen hres_1986=high_res*year_dum7
gen hres_1987=high_res*year_dum8
gen hres_1988=high_res*year_dum9
gen hres_1989=high_res*year_dum10

gen hresc_1993=mean_ptresdr_council*year_dum14
gen hresc_1994=mean_ptresdr_council*year_dum15
gen hresc_1995=mean_ptresdr_council*year_dum16
gen hresc_1996=mean_ptresdr_council*year_dum17
gen hresc_1997=mean_ptresdr_council*year_dum18
gen hresc_1998=mean_ptresdr_council*year_dum19
gen hresc_1999=mean_ptresdr_council*year_dum20
gen hresc_2000=mean_ptresdr_council*year_dum21
gen hresc_2001=mean_ptresdr_council*year_dum22
gen hresc_2002=mean_ptresdr_council*year_dum23
gen hresc_2003=mean_ptresdr_council*year_dum24
gen hresc_2004=mean_ptresdr_council*year_dum25
gen hresc_2005=mean_ptresdr_council*year_dum26
gen hresc_2006=mean_ptresdr_council*year_dum27
gen hresc_2007=mean_ptresdr_council*year_dum28
gen hresc_2008=mean_ptresdr_council*year_dum29
gen hresc_2009=mean_ptresdr_council*year_dum30

gen hrescshort_1993=mean_ptresdrshort_council*year_dum14
gen hrescshort_1994=mean_ptresdrshort_council*year_dum15
gen hrescshort_1995=mean_ptresdrshort_council*year_dum16
gen hrescshort_1996=mean_ptresdrshort_council*year_dum17
gen hrescshort_1997=mean_ptresdrshort_council*year_dum18
gen hrescshort_1998=mean_ptresdrshort_council*year_dum19
gen hrescshort_1999=mean_ptresdrshort_council*year_dum20
gen hrescshort_2000=mean_ptresdrshort_council*year_dum21
gen hrescshort_2001=mean_ptresdrshort_council*year_dum22
gen hrescshort_2002=mean_ptresdrshort_council*year_dum23
gen hrescshort_2003=mean_ptresdrshort_council*year_dum24
gen hrescshort_2004=mean_ptresdrshort_council*year_dum25
gen hrescshort_2005=mean_ptresdrshort_council*year_dum26
gen hrescshort_2006=mean_ptresdrshort_council*year_dum27
gen hrescshort_2007=mean_ptresdrshort_council*year_dum28
gen hrescshort_2008=mean_ptresdrshort_council*year_dum29
gen hrescshort_2009=mean_ptresdrshort_council*year_dum30

gen hresc_1980=mean_ptresdr_council*year_dum1
gen hresc_1981=mean_ptresdr_council*year_dum2
gen hresc_1982=mean_ptresdr_council*year_dum3
gen hresc_1983=mean_ptresdr_council*year_dum4
gen hresc_1984=mean_ptresdr_council*year_dum5
gen hresc_1985=mean_ptresdr_council*year_dum6
gen hresc_1986=mean_ptresdr_council*year_dum7
gen hresc_1987=mean_ptresdr_council*year_dum8
gen hresc_1988=mean_ptresdr_council*year_dum9
gen hresc_1989=mean_ptresdr_council*year_dum10

gen hevas_1993=high_evasion*year_dum14
gen hevas_1994=high_evasion*year_dum15
gen hevas_1995=high_evasion*year_dum16
gen hevas_1996=high_evasion*year_dum17
gen hevas_1997=high_evasion*year_dum18
gen hevas_1998=high_evasion*year_dum19
gen hevas_1999=high_evasion*year_dum20
gen hevas_2000=high_evasion*year_dum21
gen hevas_2001=high_evasion*year_dum22
gen hevas_2002=high_evasion*year_dum23
gen hevas_2003=high_evasion*year_dum24
gen hevas_2004=high_evasion*year_dum25
gen hevas_2005=high_evasion*year_dum26
gen hevas_2006=high_evasion*year_dum27
gen hevas_2007=high_evasion*year_dum28
gen hevas_2008=high_evasion*year_dum29
gen hevas_2009=high_evasion*year_dum30

gen hevas_1980=high_evasion*year_dum1
gen hevas_1981=high_evasion*year_dum2
gen hevas_1982=high_evasion*year_dum3
gen hevas_1983=high_evasion*year_dum4
gen hevas_1984=high_evasion*year_dum5
gen hevas_1985=high_evasion*year_dum6
gen hevas_1986=high_evasion*year_dum7
gen hevas_1987=high_evasion*year_dum8
gen hevas_1988=high_evasion*year_dum9
gen hevas_1989=high_evasion*year_dum10 

gen hevasc_1993=mean_ptdrevasion_council*year_dum14
gen hevasc_1994=mean_ptdrevasion_council*year_dum15
gen hevasc_1995=mean_ptdrevasion_council*year_dum16
gen hevasc_1996=mean_ptdrevasion_council*year_dum17
gen hevasc_1997=mean_ptdrevasion_council*year_dum18
gen hevasc_1998=mean_ptdrevasion_council*year_dum19
gen hevasc_1999=mean_ptdrevasion_council*year_dum20
gen hevasc_2000=mean_ptdrevasion_council*year_dum21
gen hevasc_2001=mean_ptdrevasion_council*year_dum22
gen hevasc_2002=mean_ptdrevasion_council*year_dum23
gen hevasc_2003=mean_ptdrevasion_council*year_dum24
gen hevasc_2004=mean_ptdrevasion_council*year_dum25
gen hevasc_2005=mean_ptdrevasion_council*year_dum26
gen hevasc_2006=mean_ptdrevasion_council*year_dum27
gen hevasc_2007=mean_ptdrevasion_council*year_dum28
gen hevasc_2008=mean_ptdrevasion_council*year_dum29
gen hevasc_2009=mean_ptdrevasion_council*year_dum30

gen hevascshort_1993=mean_ptdrevasionshort_council*year_dum14
gen hevascshort_1994=mean_ptdrevasionshort_council*year_dum15
gen hevascshort_1995=mean_ptdrevasionshort_council*year_dum16
gen hevascshort_1996=mean_ptdrevasionshort_council*year_dum17
gen hevascshort_1997=mean_ptdrevasionshort_council*year_dum18
gen hevascshort_1998=mean_ptdrevasionshort_council*year_dum19
gen hevascshort_1999=mean_ptdrevasionshort_council*year_dum20
gen hevascshort_2000=mean_ptdrevasionshort_council*year_dum21
gen hevascshort_2001=mean_ptdrevasionshort_council*year_dum22
gen hevascshort_2002=mean_ptdrevasionshort_council*year_dum23
gen hevascshort_2003=mean_ptdrevasionshort_council*year_dum24
gen hevascshort_2004=mean_ptdrevasionshort_council*year_dum25
gen hevascshort_2005=mean_ptdrevasionshort_council*year_dum26
gen hevascshort_2006=mean_ptdrevasionshort_council*year_dum27
gen hevascshort_2007=mean_ptdrevasionshort_council*year_dum28
gen hevascshort_2008=mean_ptdrevasionshort_council*year_dum29
gen hevascshort_2009=mean_ptdrevasionshort_council*year_dum30

gen hevasc_1980=mean_ptdrevasion_council*year_dum1
gen hevasc_1981=mean_ptdrevasion_council*year_dum2
gen hevasc_1982=mean_ptdrevasion_council*year_dum3
gen hevasc_1983=mean_ptdrevasion_council*year_dum4
gen hevasc_1984=mean_ptdrevasion_council*year_dum5
gen hevasc_1985=mean_ptdrevasion_council*year_dum6
gen hevasc_1986=mean_ptdrevasion_council*year_dum7
gen hevasc_1987=mean_ptdrevasion_council*year_dum8
gen hevasc_1988=mean_ptdrevasion_council*year_dum9
gen hevasc_1989=mean_ptdrevasion_council*year_dum10 

********************************** Figure 3A & 3B ******************************
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hres_199* hres_2000 hres_2001 hres_2002 hres_2003 hres_2004 hres_2005 hres_2006 hres_2007 hres_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hres_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hres_`i']+1.95*_se[hres_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hres_`i']-1.95*_se[hres_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)),  ylabel(0(1)3) xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("Percent") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\Fig3A.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears i.year_dum2*i.regions i.year_dum3*i.regions i.year_dum4*i.regions i.year_dum5*i.regions i.year_dum6*i.regions i.year_dum7*i.regions i.year_dum8*i.regions i.year_dum9*i.regions i.year_dum10*i.regions hres_1981 hres_1982 hres_1983 hres_1984 hres_1985 hres_1986 hres_1987 hres_1988 hres_1989 if year<=1989, fe cluster(bor_id)
qui forval i=1981(1)1989 {
	replace b_decayct_respt=_b[hres_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hres_`i']+1.95*_se[hres_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hres_`i']-1.95*_se[hres_`i'] if year==`i'
	}
replace b_decayct_respt=0 if year==1980
replace bdecayctrespt_h=0 if year==1980
replace bdecayctrespt_l=0 if year==1980
twoway (connected b_decayct_respt year if year>=1981 & year<=1989, msymbol(circle) mcolor(dkgreen) msize(medlarge) lcolor(dkgreen)) (connected bdecayctrespt_l year if year>=1981 & year<=1989, msymbol(none) lcolor(dkgreen) lpattern(shortdash))  (connected bdecayctrespt_h year if year>=1981 & year<=1989, msymbol(none) lcolor(dkgreen) lpattern(shortdash)), xlabel(1981(2)1989) legend(order(1 2) label(1 "1(High Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("Percent") xtitle("Year") graphregion(color(white))
graph save "Figures and Tables\Fig3B.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l
********************************************************************************

******************************* Figure 4 ***************************************
gen enforcement=ct_summonses/ct_reminders 
sort bor_id year
xtset bor_id year
bysort bor_id: ipolate enforcement year if ct_taxliability_dwelling!=., gen(enforc_ipolate)

** Panel A **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_taxliability_dwelling i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hres_199* hres_2000 hres_2001 hres_2002 hres_2003 hres_2004 hres_2005 hres_2006 hres_2007 hres_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hres_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hres_`i']+1.95*_se[hres_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hres_`i']-1.95*_se[hres_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)),  xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("Average tax liability per dwelling") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\Fig4A.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel B **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg enforc_ipolate i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hres_199* hres_2000 hres_2001 hres_2002 hres_2003 hres_2004 hres_2005 hres_2006 hres_2007 hres_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hres_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hres_`i']+1.95*_se[hres_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hres_`i']-1.95*_se[hres_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)),  xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("[# Court Summons]/[# Reminders]") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\Fig4B.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel C **
egen maxconctl=max(conctl) if year>=1990 & year<=1992, by(bor_id)
egen maxlabctl=max(labctl) if year>=1990 & year<=1992, by(bor_id)
egen max_conctl=max(maxconctl), by(bor_id)
egen max_labctl=max(maxlabctl), by(bor_id)

tab maxlabctl maxconctl

gen lab_1990=max_labctl*year_dum11
gen lab_1991=max_labctl*year_dum12
gen lab_1992=max_labctl*year_dum13
gen lab_1993=max_labctl*year_dum14
gen lab_1994=max_labctl*year_dum15
gen lab_1995=max_labctl*year_dum16
gen lab_1996=max_labctl*year_dum17
gen lab_1997=max_labctl*year_dum18
gen lab_1998=max_labctl*year_dum19
gen lab_1999=max_labctl*year_dum20
gen lab_2000=max_labctl*year_dum21
gen lab_2001=max_labctl*year_dum22
gen lab_2002=max_labctl*year_dum23
gen lab_2003=max_labctl*year_dum24
gen lab_2004=max_labctl*year_dum25
gen lab_2005=max_labctl*year_dum26
gen lab_2006=max_labctl*year_dum27
gen lab_2007=max_labctl*year_dum28
gen lab_2008=max_labctl*year_dum29

gen lab_1981=max_labctl*year_dum2
gen lab_1982=max_labctl*year_dum3
gen lab_1983=max_labctl*year_dum4
gen lab_1984=max_labctl*year_dum5
gen lab_1985=max_labctl*year_dum6
gen lab_1986=max_labctl*year_dum7
gen lab_1987=max_labctl*year_dum8
gen lab_1988=max_labctl*year_dum9
gen lab_1989=max_labctl*year_dum10

gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears lab_199* lab_20* i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hres_199* hres_2000 hres_2001 hres_2002 hres_2003 hres_2004 hres_2005 hres_2006 hres_2007 hres_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hres_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hres_`i']+1.95*_se[hres_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hres_`i']-1.95*_se[hres_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)),  xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("Percent") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\Fig4C.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel D **
gen nndrevasion=[ndnr_arrearsoutstanding]/ndnr_grossdebit
replace nndrevasion=nndrevasion*100

gen hres_1990=high_res*year_dum11
gen hres_1991=high_res*year_dum12
gen hres_1992=high_res*year_dum13

gen hevas_1990=high_evasion*year_dum11
gen hevas_1991=high_evasion*year_dum12
gen hevas_1992=high_evasion*year_dum13

gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg nndrevasion i.regions*i.year_dum11 i.regions*i.year_dum12 i.regions*i.year_dum13 i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 hres_199* hres_2000 hres_2001 hres_2002 hres_2003 hres_2004 if year>=1990, fe cluster(bor_id)
qui forval i=1990(1)2004 {
	replace b_decayct_respt=_b[hres_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hres_`i']+1.95*_se[hres_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hres_`i']-1.95*_se[hres_`i'] if year==`i'
	}
replace b_decayct_respt=0 if year==2005
replace bdecayctrespt_h=0 if year==2005
replace bdecayctrespt_l=0 if year==2005
sort year
twoway (connected b_decayct_respt year if year>=1990 & year<=2005, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1990 & year<=2005, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1990 & year<=2005, msymbol(none) lcolor(red) lpattern(shortdash)),  ylabel(-3(1)3) xlabel(1990(3)2008) legend(order(1 2) label(1 "1(High Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("Percent") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\Fig4D.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

grc1leg "Figures and Tables\Fig4A.gph" "Figures and Tables\Fig4B.gph" "Figures and Tables\Fig4C.gph" "Figures and Tables\Fig4D.gph", graphregion(color(white)) graphregion(margin(small))
graph save "Figures and Tables\Fig4.gph", replace
********************************************************************************


***************************** APPENDIX *****************************************
**************** Figure A1 *****************************************************
** Panel A **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears year_dum* hres_199* hres_2000 hres_2001 hres_2002 hres_2003 hres_2004 hres_2005 hres_2006 hres_2007 hres_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hres_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hres_`i']+1.95*_se[hres_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hres_`i']-1.95*_se[hres_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-1(1)4) xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA1A.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel B **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hres_199* hres_2000 hres_2001 hres_2002 hres_2003 hres_2004 hres_2005 hres_2006 hres_2007 hres_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hres_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hres_`i']+1.95*_se[hres_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hres_`i']-1.95*_se[hres_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-1(1)4) xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA1B.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel C **
gen sharepriv_1990=sharepriv*year_dum11
gen sharepriv_1991=sharepriv*year_dum12
gen sharepriv_1992=sharepriv*year_dum13
gen sharepriv_1993=sharepriv*year_dum14
gen sharepriv_1994=sharepriv*year_dum15
gen sharepriv_1995=sharepriv*year_dum16
gen sharepriv_1996=sharepriv*year_dum17
gen sharepriv_1997=sharepriv*year_dum18
gen sharepriv_1998=sharepriv*year_dum19
gen sharepriv_1999=sharepriv*year_dum20
gen sharepriv_2000=sharepriv*year_dum21
gen sharepriv_2001=sharepriv*year_dum22
gen sharepriv_2002=sharepriv*year_dum23
gen sharepriv_2003=sharepriv*year_dum24
gen sharepriv_2004=sharepriv*year_dum25
gen sharepriv_2005=sharepriv*year_dum26
gen sharepriv_2006=sharepriv*year_dum27
gen sharepriv_2007=sharepriv*year_dum28
gen sharepriv_2008=sharepriv*year_dum29

gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears sharepriv_19* sharepriv_20* lab_199* lab_20* i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hres_199* hres_2000 hres_2001 hres_2002 hres_2003 hres_2004 hres_2005 hres_2006 hres_2007 hres_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hres_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hres_`i']+1.95*_se[hres_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hres_`i']-1.95*_se[hres_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)),  ylabel(-1(1)4) xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA1C.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

grc1leg "Figures and Tables\FigA1A.gph" "Figures and Tables\FigA1B.gph" "Figures and Tables\FigA1C.gph", graphregion(color(white)) graphregion(margin(small))
graph save "Figures and Tables\FigA1.gph", replace
********************************************************************************

**************** Figure A2 *****************************************************
** Panel A **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears year_dum* hevas_199* hevas_2000 hevas_2001 hevas_2002 hevas_2003 hevas_2004 hevas_2005 hevas_2006 hevas_2007 hevas_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hevas_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hevas_`i']+1.95*_se[hevas_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hevas_`i']-1.95*_se[hevas_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-1(1)4) xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA2A.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel B **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hevas_199* hevas_2000 hevas_2001 hevas_2002 hevas_2003 hevas_2004 hevas_2005 hevas_2006 hevas_2007 hevas_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hevas_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hevas_`i']+1.95*_se[hevas_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hevas_`i']-1.95*_se[hevas_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-1(1)4) xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA2B.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel C **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears sharepriv_19* sharepriv_20* lab_199* lab_20* i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hevas_199* hevas_2000 hevas_2001 hevas_2002 hevas_2003 hevas_2004 hevas_2005 hevas_2006 hevas_2007 hevas_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hevas_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hevas_`i']+1.95*_se[hevas_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hevas_`i']-1.95*_se[hevas_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-1(1)4) xlabel(1993(3)2008) legend(order(1 2) label(1 "1(High PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA2C.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

grc1leg "Figures and Tables\FigA2A.gph" "Figures and Tables\FigA2B.gph" "Figures and Tables\FigA2C.gph", graphregion(color(white)) graphregion(margin(small))
graph save "Figures and Tables\FigA2.gph", replace
********************************************************************************


************************* Figure A3 ********************************************
** Panel A **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears year_dum* hresc_199* hresc_2000 hresc_2001 hresc_2002 hresc_2003 hresc_2004 hresc_2005 hresc_2006 hresc_2007 hresc_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hresc_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hresc_`i']+1.95*_se[hresc_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hresc_`i']-1.95*_se[hresc_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA3A.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel B **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hresc_199* hresc_2000 hresc_2001 hresc_2002 hresc_2003 hresc_2004 hresc_2005 hresc_2006 hresc_2007 hresc_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hresc_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hresc_`i']+1.95*_se[hresc_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hresc_`i']-1.95*_se[hresc_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA3B.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel C **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears sharepriv_19* sharepriv_20* lab_199* lab_20* i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hresc_199* hresc_2000 hresc_2001 hresc_2002 hresc_2003 hresc_2004 hresc_2005 hresc_2006 hresc_2007 hresc_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hresc_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hresc_`i']+1.95*_se[hresc_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hresc_`i']-1.95*_se[hresc_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("Percent") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA3C.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

grc1leg "Figures and Tables\FigA3A.gph" "Figures and Tables\FigA3B.gph" "Figures and Tables\FigA3C.gph", graphregion(color(white)) graphregion(margin(small))
graph save "Figures and Tables\FigA3.gph", replace
********************************************************************************


**************** Figure A4 *****************************************************
** Panel A **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears year_dum* hevasc_199* hevasc_2000 hevasc_2001 hevasc_2002 hevasc_2003 hevasc_2004 hevasc_2005 hevasc_2006 hevasc_2007 hevasc_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hevasc_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hevasc_`i']+1.95*_se[hevasc_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hevasc_`i']-1.95*_se[hevasc_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA4A.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel B **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hevasc_199* hevasc_2000 hevasc_2001 hevasc_2002 hevasc_2003 hevasc_2004 hevasc_2005 hevasc_2006 hevasc_2007 hevasc_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hevasc_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hevasc_`i']+1.95*_se[hevasc_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hevasc_`i']-1.95*_se[hevasc_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4)  xlabel(1993(3)2008) legend(order(1 2) label(1 "(PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA4B.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel C **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears sharepriv_19* sharepriv_20* lab_199* lab_20* i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hevasc_199* hevasc_2000 hevasc_2001 hevasc_2002 hevasc_2003 hevasc_2004 hevasc_2005 hevasc_2006 hevasc_2007 hevasc_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hevasc_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hevasc_`i']+1.95*_se[hevasc_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hevasc_`i']-1.95*_se[hevasc_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA4C.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

grc1leg "Figures and Tables\FigA4A.gph" "Figures and Tables\FigA4B.gph" "Figures and Tables\FigA4C.gph", graphregion(color(white)) graphregion(margin(small))
graph save "Figures and Tables\FigA4.gph", replace
********************************************************************************

************************* Figure A5 ********************************************
** Panel A **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears year_dum* hrescshort_199* hrescshort_2000 hrescshort_2001 hrescshort_2002 hrescshort_2003 hrescshort_2004 hrescshort_2005 hrescshort_2006 hrescshort_2007 hrescshort_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hrescshort_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hrescshort_`i']+1.95*_se[hrescshort_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hrescshort_`i']-1.95*_se[hrescshort_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA5A.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel B **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hrescshort_199* hrescshort_2000 hrescshort_2001 hrescshort_2002 hrescshort_2003 hrescshort_2004 hrescshort_2005 hrescshort_2006 hrescshort_2007 hrescshort_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hrescshort_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hrescshort_`i']+1.95*_se[hrescshort_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hrescshort_`i']-1.95*_se[hrescshort_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA5B.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel C **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears sharepriv_19* sharepriv_20* lab_199* lab_20* i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hrescshort_199* hrescshort_2000 hrescshort_2001 hrescshort_2002 hrescshort_2003 hrescshort_2004 hrescshort_2005 hrescshort_2006 hrescshort_2007 hrescshort_2008  if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hrescshort_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hrescshort_`i']+1.95*_se[hrescshort_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hrescshort_`i']-1.95*_se[hrescshort_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(Conditional PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("Percent") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA5C.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

grc1leg "Figures and Tables\FigA5A.gph" "Figures and Tables\FigA5B.gph" "Figures and Tables\FigA5C.gph", graphregion(color(white)) graphregion(margin(small))
graph save "Figures and Tables\FigA5.gph", replace
********************************************************************************

**************** Figure A6 *****************************************************
** Panel A **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears year_dum* hevascshort_199* hevascshort_2000 hevascshort_2001 hevascshort_2002 hevascshort_2003 hevascshort_2004 hevascshort_2005 hevascshort_2006 hevascshort_2007 hevascshort_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hevascshort_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hevascshort_`i']+1.95*_se[hevascshort_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hevascshort_`i']-1.95*_se[hevascshort_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA6A.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel B **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hevascshort_199* hevascshort_2000 hevascshort_2001 hevascshort_2002 hevascshort_2003 hevascshort_2004 hevascshort_2005 hevascshort_2006 hevascshort_2007 hevascshort_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hevascshort_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hevascshort_`i']+1.95*_se[hevascshort_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hevascshort_`i']-1.95*_se[hevascshort_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4)  xlabel(1993(3)2008) legend(order(1 2) label(1 "(PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA6B.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

** Panel C **
gen b_decayct_respt=.
gen bdecayctrespt_h=.
gen bdecayctrespt_l=.

xi: xtreg ct_evasion_allyears sharepriv_19* sharepriv_20* lab_199* lab_20* i.regions*i.year_dum14 i.regions*i.year_dum15 i.regions*i.year_dum16 i.regions*i.year_dum17 i.regions*i.year_dum18 i.regions*i.year_dum19 i.regions*i.year_dum20 i.regions*i.year_dum21 i.regions*i.year_dum22 i.regions*i.year_dum23 i.regions*i.year_dum24 i.regions*i.year_dum25 i.regions*i.year_dum26 i.regions*i.year_dum27 i.regions*i.year_dum28 i.regions*i.year_dum29 hevascshort_199* hevascshort_2000 hevascshort_2001 hevascshort_2002 hevascshort_2003 hevascshort_2004 hevascshort_2005 hevascshort_2006 hevascshort_2007 hevascshort_2008 if year>=1993, fe cluster(bor_id)
qui forval i=1993(1)2008 {
	replace b_decayct_respt=_b[hevascshort_`i'] if year==`i'
	replace bdecayctrespt_h=_b[hevascshort_`i']+1.95*_se[hevascshort_`i'] if year==`i'
	replace bdecayctrespt_l=_b[hevascshort_`i']-1.95*_se[hevascshort_`i'] if year==`i'
	}
replace b_decayct_respt=. if year==2009
replace bdecayctrespt_h=. if year==2009
replace bdecayctrespt_l=. if year==2009
sort year
twoway (connected b_decayct_respt year if year>=1993 & year<=2008, msymbol(circle) mcolor(red) msize(medlarge) lcolor(red)) (connected bdecayctrespt_h year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash))  (connected bdecayctrespt_l year if year>=1993 & year<=2008, msymbol(none) lcolor(red) lpattern(shortdash)), ylabel(-.1(.1).4) xlabel(1993(3)2008) legend(order(1 2) label(1 "(PT Evasion)*1(Year)") label(2 "95% CI")) ytitle("") xtitle("Year")  graphregion(color(white))
graph save "Figures and Tables\FigA6C.gph", replace
drop b_decayct_respt bdecayctrespt_h bdecayctrespt_l

grc1leg "Figures and Tables\FigA6A.gph" "Figures and Tables\FigA6B.gph" "Figures and Tables\FigA6C.gph", graphregion(color(white)) graphregion(margin(small))
graph save "Figures and Tables\FigA6.gph", replace
********************************************************************************
